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1. INTRODUCTION 

According to the system modeling, it is possible to introduce piecewise affine (PWA) systems, a 
particular hybrid system class. It is described by partitioning the vast input-state space into Polyhedron 
regions and assigning a PWA equation to each of these regions. Discrete-time PWA systems are a very 
effective tool for modeling most hybrid systems [1]. These systems have established themselves as a 
powerful class for identification and approximation of generic nonlinear systems by multi-linearization at 
equilibrium [2]. A practical method in designing controllers of nonlinear systems is optimal control concepts 
in constrained and non-constrained processes by linear discrete-time models in the form of state space. 
Study [3] presents a constrained discrete model predictive control strategy for an inside greenhouse 
temperature. Given what has been said about the modeling advantages of most systems based on the PWA 
class, recently, there has been strong interest in computing the optimal form-package controller for PWA 
systems. These problems became known as the constrained final time-optimal control (CFTOC) [4], [5]. The 
most important methods of analysis of this problem are multi-parameter programming, receding horizon 
control (RHC), or model predictive control (MPC). MPC is an effective way to deal with constrained control 
problems and has found many applications and advances in industry and academic research. This method 
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requires that the next position of the variables be predicted based on the current position, the controller input, 
and the process model. In other words, the sequence of control inputs that optimize the objective function is 
computed and applied to the process. This control concept is called the MPC [6]. 

Boiler turbine process control was performed with three manipulated variables and three controlled 
variables using MPC technique [7]. Optimal control sequence allows recalculation and feedback performance 
in MPC whenever a new measurement arrives; this is known as MPC. The stability of the control system and 
the fulfillment of conditions, constraints, and requirements during operation are ensured by the proper 
formulation of the objective function. RHC has the complexity and high volume of online computing related 
to optimizing and reducing system robustness due to the difference between the actual and MPC processes. 
Optimal methods based on the use of multi-parametric linear programming (MP-LP) and multi-parametric 
quadratic programming (MP-QP) were provided [8], [9]. 

The offline calculation of the optimal control rule for constrained discrete-time linear systems was 
performed using these methods. The resulting rules were made available as a PWA function on the 
polyhedrons. At present, explicit MPC techniques enable a standard method in controller design for nonlinear 
processes that are modeled in the PWA form, creating a substitute for intelligent controller design methods 
such as fuzzy logic and neural networks in high-performance applications [10]. Unfortunately, one of the 
main problems is the increasing complexity of the control rule obtained by increasing the prediction horizon 
and its effect on system performance. On the other hand, it is necessary to increase the prediction horizon for 
the system's optimal performance. For linear systems with parametric uncertainty by the Lyapunov function, 
the PWA controller is designed with low complexity [11]. Effective representation and approximation are 
provided by in-depth learning to MPC of linear time invariant (LTI) systems. Theoretically, at least neurons 
and hidden layers are considered [12]. A nonlinear robust MPC with input-dependent perturbations and states 
and uncertainty is presented [13]. The MPC algorithm with PWA control rules is presented for discrete-time 
linear systems in the presence of finite perturbation [14]. The online burden of computational of the linear 
MPC can be transferred to offline by using multi-parameter programming, which is named explicit 
MPC [15]. A flexibility algorithm is proposed to reduce the calculation volume in [15] that the designer can 
balance time and storage complexities. 

This is done by hash tables and the associated hash functions. Two modified controllers instead of 
the standard MP-QP are used [16] to reduce the complexity of the multi-parameter programming of MPC. 
The problem of reducing the complexity of explicit MPC for linear systems is considered by PWA 
employing separating functions [17]. A semi-continuous PWA model based on the optimal control method 
for the nonlinear system is proposed [18]. First, the nonlinear system is approximated by multi-linear 
subsystems, and then the subsystems are incorporated into a PWA system and formulated as an optimal 
control problem. A computational method for optimizing and controlling a two-tank system with three 
control valves is presented [19]. Easy implementation and optimizability of complex objective functions, 
with many local minimums are the main advantages of particle swarm optimization (PSO). Furthermore, 
PSO can be used to search significantly wider range of candidate solutions. The dynamics of the tank system 
are nonlinear. The linear model is considered, and the parameters are adjusted so that the difference between 
the actual system and the model is minimized by solving the optimal control problem. PSO has been used to 
solve the problem of constrained optimization [20]. Self-adaptive of particle swarm optimization (SAPSO) is 
recommended to increase PSO performance. Theoretically, the convergence of the method has been 
investigated. Considerable interest has recently been generated to use PSO in optimization and engineering 
problems [21]. A new algorithm that combines model predictive control with the PSO is proposed for 
optimal control of constrained direct current (DC-DC) power system modeled as piecewise affine [22]. A 
modern strategy based on model reference command shaping (MCRS) for an overhead crane, with dual 
pendulum component impacts. The existing MRCS calculation has been making strides with the PSO 
calculation to diminish plan complexity and guarantee that synchronous alteration can be done with the 
feedback controller [23]. The proposed hybrid algorithm incorporates social interaction and elitism 
mechanisms from PSO into manta ray foraging optimization (MRFO) strategy [24]. A control strategy based 
on the self-tuning method and synchronous reference frame with proportional integral (PI) regulator is 
proposed to achieve optimal power quality in an independent microgrid. Particle swarm optimization is used 
to adjust the parameters of the PI controller, which ensures flexible performance and superior power quality 
[25]. The tuning method used for proportional integral derivative (PID) with derivative filter controller for 
liquid slosh system by implementing PSO algorithm [26]. PSO algorithm is applied to optimize the surface 
grinding process parameter in both rough and final grinding conditions [27]. This paper is organized as 
follows; first, the CFTOC of PWA system is briefly described. Having introduced the two-tank optimal 
control in section 3, in section 4, the application of PSO for the solution of the expressed problem is 
discussed, and finally, the simulation results and conclusion are presented. 
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2. CONSTRAINED FINAL TIME-OPTIMAL CONTROL (CFTOC) PROBLEM AND SOLUTION 
We will do with the constrained PWA systems as (1) [4], [10]: 


x(t + 1) = fewa (x(t), u(t)) = Ajx(t) + Bult) + fi Ga) € D; (1) 


where t20, the domain D := Un 


system dynamics and D := Une D; denotes a polyhedral partition of the domain D. i.e., the closure of Dj is 
D, = f(X) € Ro |D%x + D;"u < D,°} and int(D,) Nint(D,) = @ Vi +j. The CFTOC problem for 
piecewise affine system (1) in the form (2a), (2b) and (2c) [4], [22]: 


2 D; of fpwa(.,.) is a non-empty compact set in R"**™" , Np<o is the number of 


Jr(x(0)) = min Jr(x(0), Ur) (2a) 
: x(t + 1) = fpwa(x(0), ut) 

subject to { ee (2b) 

Jr(x(0), Ur) = Ip(x(T)) + DEES 1H, uo) (2c) 


where Jr(.,.) is the cost function, 1(.,.) the stage cost, lr(.) the final penalty, Ur optimization variable 
described as the input sequence Ur = {u(t)}{i¢, T<oo receding horizon and yf is a compact terminal target 
set in R”x. If the solution of the CFTOC problem is not unique, u7(x(0)) = {u*(t)}}ig decides one 
realization from the set of the conceivable optimizers. 

The CFTOC problem determine a set of initial state and feasible inputs as xr c R™*(x(0) € xr), 
Ur- C R™ (u(t) E Up_,,t = 0,...,T — 1) respectively. The explicit closed form solution can be expressed 
as u*(t): Xr > Ur-t, t=0, ..., T-1. The considered system is PWA (1), and the cost is based on 1, œ norm 
[22]. 


U(x(t), u(t)) = Olly + IRu(e)llp (3a) 
U(x(T)) = ||Px(T)Ilp (3b) 
where ||. ||, with p={1,0} speak to the standard vector standard 1,0. The arrangement of (2) with 


aforementioned limitations is a time varying PWA work of the starting state x(O)€ P. 


u*(t) = Upwa(x(0), t) = Kr-tix (0) + Lr-ei (4) 
where t = 0,...,T, (PRE is the polyhedral partition of a set of feasible state x(0), Xr = TA Pi, with the 
closure of P; stated as P, = {x € R™*|P*x < P°} [2]. 

If a receding horizon control strategy is used for closed-loop, the control rule is expressed as 
time-variant PWA state-feedback of the form [4]: 


Men (x(t)) = Krix(t) + Lr; if x(t) € P; (5) 


where i=1,..., Np and for t20, u*(t) = ugy (x(t). The CFTOC problem can be presented and solved for any 
selection of P, Q, R, albeit here it is assumed that the parameters T, Q, R, P and yf are selected by the 
following assumptions [3]. To avoid additional control actions in steering states to the origin (equilibrium 
point), matrices R, Q are required to have a full column rank. 


3. OPTIMAL CONTROL OF TWO-TANK SYSTEM 

The two-tank [28] shown in Figure 1 is a basic benchmark model to investigate and analyze the 
control issues for PWA system. The tanks are filled by pump acting on tank 1, continuously manipulated 
from 0 up to a maximum flow Q1. A switching valve V12 controls the flow between the tanks. This valve is 
assumed to be either completely opened or closed (V12=0 or 1 respectively). The Vn2 manual valve controls 
the nominal outflow of the second tank. It is assumed in the simulations that the manual valves, Vni is always 
closed and that Vn2is open. The liquid levels to be controlled are denoted by hi, hz for each tank respectively. 
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Figure 1. Two-tank system schematic 


It consists of two tanks that connected to each other. We assumed that: i) the valves behavior is 
linear, ii) the initial volume of liquid in tanks is zero, and iii) the inflow of liquid to the first tank is constant 
and has the maximum value. The liquid volume of tank 1 is defined by a time varying equation as (6a) and 
(6b): 

V1(t) = V0,1 + (Q1 — Q1,2) xt (6a) 

V1(t) = Al x h1(t) (6b) 
where Vo, is initial volume of liquid in tank land Q1, Qı, 2 are inflow and outflow liquid of tank! (Q1,2 can be 
defined as inflow of liquid to tank 2), A; and h; are base area and the time varying height of liquid in tank1 
respectively. The (6) can be repeated for the tank 2 with similar definition as (7a) and (7b): 

V2(t) = V0,2 + (Q1,2 — Q2) xt (7a) 

V2(t) = A2 x h2(t) (7b) 


By combining (6), (7), 
h(t) = oa + (Q = Q12) x t) (8a) 


1 
h(t) = a, o2 + (Q1,2 — Q2) x t) (8b) 
The system is expressed as a discrete time model with a sample time (T,;=10s) by (9). 


hy (k + 1) = hy (k) +2 (Q1 (k) — Kex2Vin (Ia (k) — fa (k)) 


Ay 


hy (ke + 1) = ha (Kk) + Ë (hea2Vii (Ita (k) — ha(k)) — KwaVivaha (Ie) 


A2 


(9) 


This model can be defined as a piecewise affine system of the shape (1), with four modes, depicted as takes 
after [28]: 

— First Mode: Vř, open, hi>hy 

— Second Mode: Vř, open, hi<hy 

— Third Mode: Vý, closed, hi=hy 

— Fourth Mode: Vý, closed, hi<hy 
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0.0164 
—0.0164 


i= 


0.9542 —0.0393) , _ [0.0699 07, _ a _ 
ee 0.9670 |2 =| 0 al C, = [10] D, = [0] and f = | 


The CFTOC problem of the presented PWA system is solved by MPT [29] based on MPC for the 
prediction horizon=3, norm =1, Q = p | R= 10™ x f `] and the explicit PWA control law has 78 
polyhedral regions as shown in Figure 2 and the performance of the closed-loop system from a given initial 
condition is presented in Figure 3. 

Using PSO algorithm, the considered purposes are going to be fulfilled simultaneously: 

— The number of polyhedral of explicit MPC-based control law is minimized to reduce the complexity, 
— The liquid reaches a certain height in tanks in a short time and desirable manner. 


Controller partition with 42 regions 


Level 2 


0 0.1 0.2 0.3 0.4 0.5 0.6 
Level 1 


Figure 2. Controller partitions 
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Figure 3. Closed loop system performance 
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4. PSO ALGORITHM APPLICATION FOR SOLVING DEFINED PROBLEM 

Before long a short time later, a brief audit of PSO calculation is displayed, and after that PSO 
application is explored to fathom the defined problem. The PSO could be a worldwide exploratory 
optimization strategy, to begin with, proposed by Kennedy and Eberhart in 1995 [30]. It is made of swarm 
insights and is based on inquiring about the behavior of feathered creatures and the angle of development. 
The motion of each molecule is impacted by its best-known nearby position and is additionally coordinated 
to the best-known positions within the look space, which are upgraded by other particles to discover superior 
positions. Agreeing to Figure 4, the premise of the strategies is as follows: “Each molecule can appear with 
its current speed and position, the foremost ideal position of each individual, and the foremost ideal position 
of the environment” [22], [30]. By selecting the beginning populace Xj, Vi, the speed and position of each 
molecule around the look space alter concurring to correspondence (10) [30]: 


Vie = VÉ +c x rë x (VP) + c rt x (VE) (10a) 
Xia" = Xia + Vig (10b) 
vest = pbest% — XE, (10c) 
Vc = ghest — X% (10d) 


where in this uniformity, VS and X% independently stand for the speed of the molecule i at its k times and the 
d-dimension amount of its position; pbest% speaks to the d-dimension amount of the person i at its most 
ideal position at its k times; pbest*, is the amount d-dimension of the swarm at its most ideal position. To 
avoid a molecule from being distant from the look space, the speed of the made molecule is restricted in each 
course between -Vamax and Vamax [22]. In the event that the number of Vamax is as well tall, the arrangement is 
distant from the most excellent arrangement, something else; the ideal arrangement will be nearby. cı and c2 
appear the speed digit, which alters the length when flying to the biggest molecule within the entire gather 
and the ideal single molecule. In the event that the number of Vamax is as well little, the molecule is likely 
distant from the target field; in case the shape is as well huge, the molecule may fly all of a sudden to the 
target field or past the target field. 


Figure 4. The basis of evolutionary PSO algorithm 


Reasonable figures for cı and c2 can control the molecule flight speed and will not be a halfway ideal 
arrangement. cı is more often than not rise to c2, and they are rising to 2. rı and r2 speak to an irregular story 
and 0-1 could be an arbitrary number. As said some time recently, our modern point is utilizing PSO for 
complexity lessening of unequivocal MPC-based control law by diminishment the number of its polyhedral 
and setting the physical parameter of framework to move forward the framework execution at the same time. 
In this manner, the taking after objective work has been defined: Fitness-Function =Number of polyhedral + 
Output specifications. 
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Where yield determinations are decided as a summation of operational determinations such as 
settling time, overshoot, undershoot, steady-state deviation, time consistent, and so forward. Presently, the 
unequivocal controller gotten in past segment as a portion of PSO ought to be considered and the taking after 
unused execution list is being defined: 


Jnew = Min(Fitness — Function) = Min[number of polyhedrals of ury(x(0) 
+ Output speciication] (11a) 


rf" = URH (x(t) = Krix(t) + Ly; ifx(t) E P; (11b) 


Output = )' desired output characteristic 


It is used according to the flowchart shown in Figure 5. 


Figure 5. Flowchart of PSO algorithm application 
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The issue is illuminated as: 

— The required information such as time step, fluid stature, and calculation parameters such as the number 
of population and a number of the cycles is considered. 

— The initial population is made as follows [22]: 


X; = [height of valve1, cross section of valves, maximum height, base area of each tank] 
Xı 

eee , X2 

initial population = ; 


XN pop 


— The yield determinations of the framework are measured based on introductory parameters, at that point, 
the control law (gotten by MPT) is connected, and the number of polyhedral districts is calculated. In the 
long-run objective work is characterized as [22]: 


Fitness — Function = Number of polyhedral + output specifications. 


where we consider yield determinations can be expected as yield details £ settling time + certain stature 
of liquid. The best arrangement among the entire population is decided, and the population is overhauled 
based on (10a) and (10b). 
— For the predefined number of iterations, step 3 is done repetitively. 
The convergence conditions are explored, and at last the leading arrangement is appeared within the 
output. Then number of iterations is 15, the number of populaces is 5, cı= c2=2 is considered. The ideal 
parameters are compiled within the Table 1. 


Table 1. The optimal parameters 


Sampling time _Height of valve |__Cross-section of valves | Maximumheight Areaofeachtank NumP __ Time(S) 
10s 0.1 m 1.00E-05 0.5 0.001 m^2 10 189 


T is the required time to reach a certain height of the liquid. As presented, the number of control law 
polyhedral reduces from 42 to 10. In Figure 6, the output flow variation of the first tank is shown. These 
changes are linear. After 100 seconds, output flow increases to 10 cubic meters. In Figure 7, the liquid height 
in the first tank is shown. These changes are nonlinear and incremental. After 90 seconds, the liquid height of 
the first tank reaches 50 meters. Figure 8 shows the output flow variation of the second tank. These changes 
are linear. At the last moment, there was a change in the function of the output valve. 


Across Section Of Valves Steady State Level In First Tank 


Water Flow Rate nD 


10 


0 50 100 150 300 o 50 100 150 200 
Time (S) Time (S) 
Figure 6. The output flow variation of tank 1 Figure 7. The liquid height in tank 1 


Figure 9 shows the liquid height changes of the second tank. As you can see, this change is 
nonlinear and incremental. The height of the liquid in both tanks reached a certain value in 189 seconds. 
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Figure 10 shows the unique convergence of the objective function using the PSO algorithm. It converges 
after 4 iterations. Figure 11 shows the controller partition with 10 regions. Concurring to the comes about, it 
can be concluded that by tuning the physical parameters utilizing PSO, the number of polyhedral of 
MPC-based control law is minimized; in this way, the related complexity of arrangement to CFTOC will be 
diminished together with the change of yield details at the same time so that the fluid tallness in both tanks 
reach craved esteem in 189 seconds. 


Output Valve Of Second Tank 
10p 7 i F + 60; 


Steady State Level In First Tank 


7 F 


Water Flow Rate w) 
wn 
è 


Height Of Water (m) 
w 
ô 


Y 
i=} 
T 


10}- 


rt E 5 PE ais $ = —— 
0 50 _ 100 150 200 % 50 100 150 200 
Time (S) Time (S) 
Figure 8. The output flow variation of tank 2 Figure 9. The liquid height in tank 2 
45r : p p r r ~ Controller partition with 10 regions 


Number Of Polyhedral 
Level 2 


2 4 6 8 10 12 14 0 0.1 0.2 0.3 0.4 0.5 
Iteration Level 1 


Figure 10. Convergence diagram of objective Figure 11. Controller with minimum partitions 
function 


5. CONCLUSION 

Several analytical methods have been used for the CFTOC solution. Their main disadvantage is the 
computational complexity of the solution; therefore, the problem can be considered as NP-hard. The meta- 
heuristic algorithm is used to solve NP-hard optimization problems that have strategies to escape the local 
optimal solution and are applied to a wide range of problems. In general, the developments of meta-heuristic 
methods are done with the study and inspiration of optimization in nature, such as particle swarm 
optimization. By using PSO and an appropriate definition of the objective function, the complexity of the 
MPC-based solution of CFTOC was reduced and the system performance was improved simultaneously. The 
most massive advantage of the recommended method is that if the mentioned purposes were not in one 
direction, we can define a multi-objective function to fulfill aims. According to the simulation results, it is 
demonstrated that the number of polyhedral and the dependent complexity of CFTOC solution are reduced, 
the system performance such as reaching the liquid height at a certain time is desirable and the obtained 
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steady-state error reaches zero. The number of control law polyhedral reduces from 42 to 10. The liquid 
height in both tanks reaches desired value in 189 seconds. 
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